A&A manuscript no. 

(will be inserted by hand later ] 



Your thesaurus codes are: 

08 (08.01.03; 08.06.03; 08.12.02; 08.09.02) 



ASTRONOMY 

AND 
ASTROPHYSICS 



On the Interpretation of the Optical Spectra of L-type 
Dwarfs* 



o 
o 
o 

(N 
G 

in 



> 
o 
\o 
o 

o 
o 
o 



Oh' 
i 1 

O 1 



X 



Ya. Pavlenko 1 , M. R. Zapatero Osorio 2 , and R. Rebolo 2 3 

1 Main Astronomical Observatory, Golosiiv, 03680 Kyiv-127, Ukraine 

2 Instituto de Astrofi'sica de Canarias. C/. Via Lactea s/n, E-38200 La Laguna, Tenerife. Spain 

3 Consejo Superior de Investigaciones Cientfficas, Spain 

Received ; accepted 



Abstract. We present synthetic optical spectra in the red 
and far-red (640-930 nm) of a sample of field L dwarfs 
suitably selected to cover this new spectral class, and the 
brown dwarf GL229B. We h ave used the rece nt "d usty" 
atmospheres by Tsuji ( 2000 ) and by Al lard ( 1999 ), and 



a synthesis code (Pavlenko et al. 199E ) working under 
LTE conditions which considers the chemical equilibrium 
of more than 100 molecular species and the detailed opac- 
ities for the most relevant bands. Our computations show 
that the alkali elements Li, Na, K, Rb, and Cs govern 
the optical spectra of the objects in our sample, with Na 
and K contributing significantly to block the optical emer- 
gent radiation. Molecular absorption bands of oxides (TiO 
and VO) and hydrides (CrH, FeH and CaH) also domi- 
nate at these wavelengths in the early L-types showing a 
strength that progressively decreases for later types. We 
find that the densities of these molecules in the atmo- 
spheres of our objects are considerably smaller by larger 
factors than those predicted by chemical equilibrium con- 
siderations. This is consistent with Ti and V atoms being 
depleted into grains of dust. 

In order to reproduce the overall shape of the optical 
spectra of our observations an additional opacity is re- 
quired to be implemented in the computations. We have 
modelled it with a simple law of the form a (v/v ) N , with 
N = 4, and found that this provides a sufficiently good fit 
to the data. This additional opacity could be due to molec- 
ular/dust absorption or to dust scattering. We remark 
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that the equivalent widths and intensities of the alkali 
lines are highly affected by this opacity. In particular, the 
lithium resonance line at 670.8 nm, which is widely used 
as a substellarity discriminator, is more affected by the 
additional opacity than by the natural depletion of neu- 
tral lithium atoms into molecular species. Our theoretical 
spectra displays a rather strong resonance feature even at 
very cool effective temperatures (~1000K); depending on 
the effective temperature and on the amount of dust in the 
atmospheres of very cool dwarfs, it might be possible to 
achieve the detection of lithium even at temperatures this 
cool. Changes in the physical conditions governing dust 
formation in L-type objects will cause variability of the 
alkali lines, particularly of the shorter wavelength lines. 



Key words: Stars: atmospheres; fundamental parame- 
ters; low-mass, brown dwarfs; individual: BRI 0021-0214, 
G1229B, Kelul, DenisP J1228-1547, DenisP J0205-1159 



1. Introduction 

The optical spectra of the recently discovered very cool 
dwarfs present new challenges to theoretical interpreta- 
tion. Their spectral characteristics are drastically differ- 
ent from those of the well known M-dwarfs and this has 
prompted the use of a new spectral classification, the so- 
called L-dwarfs (Martin et al. |l997a| |1999| ; Kirkpatrick 
et al. 1999a). In principle, the study of the optical spec- 
tral energy distribution may allow a better understand- 
ing of their physical properties, effective temperatures, 
gravities, atmospheric composition, etc. The main molec- 
ular absorbers at optical wavelengths in early- to mid-M 
dwarfs are significantly depleted at the lower tempera- 
tures present in L dwarfs because of the incorporation of 
their atoms into dust grains. This process should start 



in the latest M-dwarf atmospheres (Lunine et al. 1989 



Tsuji, Ohnaka & Aoki 1996a; Tsuji et al. 1996b; Jones 



Tsuji 1997; Allard et al. 1997), considerably reducing 



the strength of TiO, VO and other molecular bands, and 
producing significant changes in the overall properties of 
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Table 1. Log of spectroscopic observations 



Object 




Spec. Type 


Telescope 


Dispersion 


AA 


Date 


Exposure 










(A/pix) 


H 


(1997) 


(s) 


Kelul 




L2 (L2) 


WHT 


2.9 


640-930 


Jun 16 


2x1200 


Denis-P J 1228- 


-1547 


L4.5 (L5) 


WHT 


2.9 


640-930 


Jun 14 


2x1200 


Denis-P J0205- 


4159 


L5 (L7) 


WHT 


2.9 


640-930 


Aug 23 


2x1800 








Keckll 


2.5 


640-890 


Nov 2 


1x1600 



NOTES. Spectral types are given in Martin et al. (1999). Those spectral types in brackets come from Kirkpatrick et al. (1999a 



the optical spectra. Naturally, the appearance of dust will 
modify the temperature structure of the atmosphere, sig- 
nificantly affecting the formation of the emerging spec- 



Imaging Spectrograph (LRIS, Oke et al. 1995| ) with the 



tram (Allard et al. |1997| ). 

In this paper we follow a semi-empirical approach to 
understand the relevance of different processes on the 
resulting spectral energy distributions in the optical for 
L dwarfs and G1229B. We have obtained far-red opti- 
cal spectra of several of these cool dwarfs and compared 
them to synt hetic spectra gener ated using the latest mod- 
els by Tsuji J2000[ ) and Allard ( |l999| ). These models have 
been successfully used in the interpretation of the near- 
infrared spectra of these objects (Tsuji, Ohnaka & Aoki 



199S ; Kirkpatrick et al. 1999b ). We have taken into ac 



count the depletion of some relevant molecules associated 
with the formation of dust and we have investigated the 
effects of dust scattering and/or absorption on the for- 
mation of the optical spectra. Remarkably strong alkali 
lines are present in the spectra and dominate its shape 
in the 600-900 nm region, providing major constraints to 
the theoretical modelling. We present the observations in 
section 2, models and synthetic spectra in section 3, and 
the role of alkalis, molecular bands and dust scattering 
and/or absorption is considered in section 4. In section 5 
we discuss the implications of our study on the determi- 
nation of effective temperatures and gravities as well as 
the formation of the lithium resonance line which is a key 
discriminator of substellar nature in these objects. Finally, 
conclusions are presented in Section 6. 



2. Observations and data reduction 

We have collected intermediate-resolution spectra in 
the 640-930 nm range for Kelul (Ruiz, Leggett & Al- 
lard [1997D , Den is-P J 1228-1547 and Denis-P J0205-1159 
(Delfosse et al. 1997) using the 4.2m William Herschel 
Telescope (WHT; Observatorio del Roque de los Mucha- 
chos, La Palma). The coolest L dwarf in our sample, Denis- 
P J0205-1159, has also been observed with the Keckll tele- 
scope (Mauna Kea Observatory, Hawaii). Table [l] sum- 
marizes the log of the observations. The instrumentation 
used was the ISIS double-arm spectrograph at the WHT 
(only the red arm) with the grating R158R and a TEK 
(1024x1024 pix 2 ) CCD detector, and the Low Resolution 



600 grove mm" 1 grating and the TEK (1024x 1024 pix 2 
CCD detector at the Keckll telescope. The nominal dis- 
persion and the wavelength coverage provided by each in- 
strumental setup were similar and are listed in Table [l]. 
Slit projections were typically 2-3 pix giving spectral res- 
olutions of 6-8 A. Spectra were reduced by a standard 
procedure using IRAffJ which included debiasing, flat- 
fielding, optimal extraction, and wavelength calibration 
using the sky lines appearing in each individual spectrum 
(Osterbrock et al. 1996 ). Finally, the observed WHT spec- 
tra were corrected from instrumental response making use 
of the spectrophotometric standard stars BD +26° 2606 
and HD 19445 which have absolute flux data available in 
the IRAF environment. No flux standards were observed 
at the Keckll telescope, and thus the instrumental signa- 
ture of the Keckll spectrum of Denis-P J0205-1 159 was 
removed by matching it with the WHT spectrum. Our 
observations are displayed in Fig. |l| together with spec- 
tra of the other two objects in our sample: BRI 0021-0214 
(M9.5, Irwin, McMahon, & Reid |1991[ ), and G1229B (a 
brown dwarf with methane absorptions, Nakajima et al. 
|1995 ; Oppenheimer et al. 1995). Our sample covers a wide 
range of spectral types, from the transition objects be- 
tween M and L types down to the latest types. 



3. Model atmospheres, spectral synthesis code, 
chemical equilibrium and opacities 

We carried out the computations using the LTE spectral 
synthesis program WITA5, which is a modified version of 
the program used by Pavlenko et al. (1995) for the study 
of the formation of lithium lines in cool dwarfs. The mod- 
ifications were aimed to incorporate "dusty effects" which 
can affect the chemical equilibrium and radiative trans- 
fer processes in very cool atmospheres. We have used the 
set of Tsuji's ( |2000[ ) "dusty" (C-type) LTE model atmo- 
spheres. These models were computed for the case of segre- 
gation of dust-gas phases, i.e. for conditions r^ust > r cr it, 
where Tdust is the size of dust particles and r cr i t is critical 



1 IRAF is distributed by National Optical Astronomy Obser- 
vatories, whcih is operated by the Association of Universities 
for Research in Astronomy, Inc., under contract with the Na- 
tional Science Foundation. 
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Fig. 1. Far-red optical spectra for our sample. New observations are those of Kelu 1, Denis-P J1228-1547 and Denis- 
P J0205-1159 (Keckll spectrum). Data for BRI 0021-0214 and G1229B have been collected from Martin, Rebolo & 
Zapatero Osorio (1996) and from Schultz et al. (1998), respectively. Spectral types for the L-dwarfs are given following 
the classification of Martin et al. (1999). Spectra in the left panel have been shifted by 0.8 units for clarity. Identification 
of some atomic and molecular features is provided in the top. 



size corresponding to the gas-dust detailed equilibrium 
state. In a previous study (Pavlenko, Zapatero Osorio & 



Rebolo |2000D we had used Tsuji's ( |2000D B-type models 
which are computed for the case of Tdust = r C rU- In this 
paper we have also used a grid of the NextGen " dusty " 
model atmospheres computed recently by Allard ( |l999|) . 
The temperature-pressure stratification of Allard's mod- 
els lie betwe en those of the C-type and B-type models of 
Tsuji ( [2000 ) as it can be seen in Fig. [| 

Chemical equilibrium was computed for the mix of 
«100 molecular species. Alkali-contained species forma- 
tion processes were considered in detail because of the 
important role of the neutral alkali atoms in the forma- 
tion of the spectra. Constants for chemical equil ibrium 
computations were taken mainly from Tsuji ( 1973 ). 



In the high pressure conditions of the atmospheres of 
L-dwarfs some molecules can be oversaturated (Tsuji et 
al. 1996a); in this case such molecules should undergo con- 



densation. To take into account this effect, we reduced the 
abundances of those molecular species down to the equi- 



librium values (Pavlenko 1998). The constants for compu- 
tations of saturation densities were taken from Gurwitz et 
al. 



(1982) 



We used the set of continuum opacity sources listed 
in Table |^ where we also give the original sources for the 
opacity computation codes. That opacity grid allows us to 
obtain resonable fits to a variety of stars (see e.g. Martin et 
al. 1997t; Israelian, Garcia Lopez & Rebolo 1998; Yakov- 
ina & Pavlenko 1998). Opacities due to molecular band ab- 
sorption were treated using the Just Overlapping Line Ap- 
proximation (JOLA). Synthetic spectra of late M-dwarfs 
using both the continuum opacities listed in Table ^| and 
the JOLA treatment for molecular b and a bsorptions have 
been already discussed in Pavlenko ( |l997|) . 



The alkali line data were 
database (Piskunov et al 



taken from the VALD 
1995|, see their Table 3). At 



the low temperatures of our objects we deal with satu- 
rated absorption lines of alkalis. Their profiles are pressure 
broadened. At every depth in the model atmosphere the 
profile of the aborption lines is described by a Voigt func- 
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Table 2. List of the opacity sources used in our computations 



Opacity source 



Subprogram Author 



Bound-free absorption of H HOP 

absorption H2PLOP 

ET ion absorption HMINOP 

Rayleigh scattering by H atoms HRAYOP 

He I bound-free absorption HEIOP 

Hell bound- free absorption HE20P 

He" absorption HEMIOP 

Rayleigh scattering by He HERAOP 
Bound-free absorption ofMgl, All, Si I, Fel, OH and CH COOLOP 

Bound-free absorption of Ni, Oi, Mgn, Sin, Cai LUKEOP 

Bound- free absorption of C II + N II + O II HOTOP 

Electron Thompson scattering by e~ ELECOP 

Rayleigh scattering by H 2 H2RAOP 

Absorption due to Hydrogen lines, HLINOP 

H2 quasimolecule absorption QUASIH 

CT ion absorption CMINUSOP 

Bound- free of C 1 ClOPAC 



Kurucz ( 1993 ) 
Kurucz ( 1993 ) 
Kurucz ( 1993 ) 
Kurucz ( 1993 ) 
Kurucz ( 1993 ) 
Kurucz ( 1993 ) 
Kurucz ( 1993 ) 
Kurucz ( 1993 ) 
Kurucz ( 1993 ) 
Kurucz ( 1993 ) 
Kurucz ( 1993 ) 
Kurucz ( 1993 ) 
Kurucz (1993) 
Kurucz (1993) 
Pavlenko ( h999a 
Pavlenko ( 1999a 
Pavlenko ( 1999a 
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Table 3. List of atomic lines used in our computations 



Pressure (din cm 2 ) 

Fig. 2. Temperature structures for the model atmo- 
spheres of T e g = 1400K and log g = 5.0 given by Tsuji 
( 200C , C-typ e, dashed line; B-type, dotted line) and by 
Allard ( |l999i full line). 



tion H(a,v), where damping constants a were computed 
as in Pavlenko et al. ( 1995 ). 



3.1. Molecular opacities 

Detail molecular opacities have been computed for all the 



molecules listed in Table 2 of Pavlenko et al. (1995) as 
well as for CrH and CaH. To compute the opacity due to 
absorption of VO and TiO bands we followed the scheme 





X vac (nm) 


gf 


E" (eV) 


Li 1 


610.5164 


0.126 E+01 


1.840 


Li 1 


610.5275 


0.229 E+01 


1.850 


Li 1 


610.5290 


0.257 E+00 


1.850 


Li 1 


670.9551 


9.790 E-01 


0.000 


Li 1 


670.9701 


4.900 E-01 


0.000 


Li 1 


812.8408 


0.216 E+00 


1.850 


Li 1 


812.8629 


0.432 E+00 


1.850 


Nal 


589.1518 


1.288 E+00 


0.000 


Nal 


589.7489 


6.457E-01 


0.000 


Nal 


819.6986 


3.390 E-01 


2.105 


Nal 


819.7016 


3.090 E+00 


2.105 


Nal 


818.5443 


1.690 E+00 


2.100 


Ki 


766.6961 


1.340 E+00 


0.000 


Ki 


770.1031 


6.760 E-01 


0.000 


Rbi 


780.2348 


1.370 E+00 


0.000 


Rbi 


794.9729 


6.800 E-01 


0.000 


Csi 


894.5900 


7.800 E-01 


0.000 


Csi 


852.3285 


1.620 E+00 


0.000 



presented in Pavlenko et al. (1995) and Pavlenko (1997). 
However, for the £? 4 II( r ) — X 4 S _ band system of VO we 
used a more complete matrix of Franc-Condon Factors 



computed by the FRANK program (Cymbal 1977) with 
account of rotational- vibrational interaction in the Morse- 
Pekeris approximation modified by Schumakcr ( |1969 



sec 



Pavlenko 1999b for more details). Futhermore for this 
band system we used oscillator strength f e from Allard 



Hauschildt (1995). Fo r the e band system of TiO the 



parameters of Schwenke ( 1995 ) were used. 

The CrH molecular band opacity was also consid- 
ered, the data required for the computations of the elec- 
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tronic transition A 6 £( + )-X 6 £( + ) being taken from Hu- 
ber & Herzberg (1979). Franc-Condon factors were com- 
puted by the FRANK program (Cymbal |1977| ). For this 
band system we used f e = 0.001, determined by Pavlenko 
( |!999bp . Molecular bands of the A 6 £(+)-X 6 £(+) of CrH 
lie in the wide wavelength region 500-1200 nm, and a head 
of the strong (0,1) band lies near the core of the K I atomic 
line. A head of the (0,0) band of A 6 S]M-X 6 S]M system 
of CrH at 860 nm is blended with the heads of (1,0), (3,2), 
(2,1) bands of the B 4 n (r) -X 4 £- system of VO. For the 
spectral region 650-710 nm we took into account the ab- 
sorption due to the B 2 T,-X 2 T, band system of CaH for 
which we adopted / e = 0.05. Frank-Condon factors were 
computed with the FRANK program. Although FeH is 
an important absorber at wavelengths around 870 nm and 
990 nm for the coolest dwarfs, it is not yet included in 
our calculations due to the lack of appropriate laboratory 
data. In the wavelength region presented in this paper, the 
blue band of FeH is blended with CrH and VO. 

4. Analysis and interpretation 

4--1- Atomic features 

Our spectral synthesis for model atmospheres in the tem- 
perature range 2200-1000 K confirms the relevant role of 
atomic features due to the Na I and K I resonance doublets 
in the spectral range 640-930 nm. The strength of these 
doublets increases dramatically when decreasing the ef- 
fective temperature (T cf f). In Fig. || we display synthetic 
spectra for different T e g and gravity values. These com- 
putations do not include any molecular or grain opacity 
in order to show how alkali absorptions change with these 
parameters. In addition, we can see in the figure how the 
overall shape of the coolest L-dwarf spectra is governed by 
the resonance absorptions of Nal and Ki. Even the less 
abundant alkali (i.e. Li, Rb, Cs) produce lines of remark- 
able strength. Note the increase in intensity of K I and Na I 
lines with decreasing T g and with increasing atmospheric 
gravity. Although this different behaviour can, in princi- 
ple, make it difficult to disentangle these parameters from 
the optical spectra, simple physical considerations give an 
upper limit to gravity of \ogg = 5.5 for brown dwarfs with 
lithium and therefore the large broadening of the K I lines 
seen in some of our objects cannot be attributed to higher 
aphysical values of gravity. 

Our computations provide a qualitative explanation 
of the far-red spectral energy distributions presented in 
Fig. |l|. The equivalent widths (EWs) of the Ki and 
Nal resonance doublets may reach several thousand 
Angstroms, becoming the strongest lines so far seen in 
the spectra of astronomical objects. This is mainly caused 
by the high pressure broadening in the atmospheres of the 
coolest dwarfs, where damping constants of the absorption 
lines vary from 0.001 in the uppermost layer of the atmo- 
sphere to 2-4 in the deepest regions. The subordinate lines 
of Nal at 819.5 nm, clearly seen in all the L-dwarfs in our 
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Fig. 3. Theoretical optical spectra for the alkali ele- 
ments computed with different values of T e g and grav- 
ity. The upper panel shows the dependence with temper- 
ature (computations performed using Tsuji's 2000 C-type 
models and \ogg = 5.0. Spectra are normalized at 880 nm). 
The lower panel displa ys th e dependence with gravity at 
T eff = 1600K (Allard's |i~999| "dusty" models. Spectra are 
normalized at 860 nm). Identifications of the atomic lines 
are given in the top. 



sample, become weaker as T e g decreases. These computa- 
tions also show that the subordinate Lil line at 812.6 nm 
may be detected in early/mid L-dwarfs with equivalent 
widths not exceeding lA, and that the triplet at 610.3 nm 
appears completely embedded by the wings of the K I and 
Nal resonance lines. 
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Fig. 4. The observed spectrum of Kelu 1 normalized at 
~840 nm is shown by the full line and the predicted spec- 
tra (Tsuji's C-type model atmosphere for T e ff = 2000K 
and log<? = 5) are shown by the dash-dotted line at a res- 
olution of 10A. Upper panel displays the theoretical spec- 
trum computed considering only the natural depletion of 
molecules as a result of chemical equilibrium. Lower panel 
depicts the same computations taking into account an 
"extra" depletion of TiO (# = 0.01), VO (# = 0.08), CaH 
(# = 0.5) and CrH (# = 0.2). 



Fig. 5. The observed spectrum of DenisP J1228-1547 nor- 
malized at ~860nm is shown by the full line and the 
predicted spectra (Tsuji's C-type model atmosphere for 
T e g = 1600 K and log g = 5) are shown by the dash-dotted 
line at a resolution of 10 A. Upper panel displays the the- 
oretical spectrum computed considering only the natural 
depletion of molecules as a result of chemical equilibrium. 
Lower panel depicts the same computations taking into ac- 
count an "extra" depletion of TiO (# = 0.0), VO (# = 0.0), 
CaH (# = 0.5) and CrH (# = 0.04). 



4-2. Molecular features 

In the optical spectra of very cool dwarfs one expects the 
presence of bands of VO, TiO and indeed, our syntheti- 
cal spectra show that these bands play an important role 
in the cases of BRI 0021-0214 and Kelu 1. In Fig. g (up- 
per panel) we compare the observed spectrum of Kelu 1 
with a synthetic spectrum obtained using the Tsuji C- 
type model for T o ff = 2000K and log<? = 5. The natural 
depletion of molecules resulting from the chemical equi- 



librium is not sufficient to get a reasonable fit to the data. 
The discrepancies with respect to the observed spectrum 
can be notably reduced (Fig. lower panel) if we impose 
a depletion of CaH, CrH, TiO and VO molecules which 
should account for the condensation of Ca, Cr, Ti and 
V atoms into dust grains. Since we do not have an ap- 
propriate theoretical description of the processes of dust 
formation at present, we use the simple approach of modi- 
fying the chemical equilibrium for all molecules. We imple- 
mented this "extra" depletion simply by introducing a fac- 
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Fig. 6. Kelul's observed spectrum (full line) compared 
t o the best fits (dash-dotted line) obtained using Tsuji's 
fl2000p C-type models (uppe r panel, T off = 2000K and 
log<? = 5), and Allard's ( 1999| ) dusty models (lower panel, 
T e ff = 1800K and logg = 5). Both predicted spectra with 
a resolution of 10A have been computed considering the 
R depletion factors given in Fig. ^ and the AdO law de- 
scribed in the text with a Q =0.03. 



Fig. 7. DenisP J1228-1547 and DenisP J0205-1159's ob- 
served data (full line) compared to the best fits (dash- 
dotted line) obtained using Tsuji's fl2000| ) C-type mod- 
els (upper panel: T e ff = 1600K and logg = 5; lower panel: 
T e ff = 1200K and logy = 5). Both predicted spectra with 
a resolution of 10A have been computed considering the 
R depletion factors given in Fig. || and the AdO law de- 
scribed in the text with a a =0.006. 



tor R which describes the reduction of molecular densities 
of the relevant species over the whole atmosphere. From 
the comparison between observed and computed spectra 
we find that the reduction factor for TiO ranges from 
(complete depletion) to 1 (i.e. non depletion). In Fig. || 
we can see the comparison of similar synthetic spectra 
(Tsuji's C-type model atmosphere for T o ff = 1600K and 
log# = 5) with the spectrum of DenisP J1228-1547. Total 
or almost total depletion of Ti and V into the dust grains is 
required to explain the spectrum of the mid- type L-dwarfs 
DenisP J1228-1547 and DenisP J0205-1159. 



4-3. The need for additional opacity (AdO) 

Our first attempts to model the spectra including atomic 
and molecular features showed in Figs. ^ and ^| were only 
modestly successful. Although we could reproduce reason- 
ably well the red wing of the K I line, the theoretical fluxes 
in the blue part of our synthetic spectra (640-750 nm) are 
too large. In fact, we cannot fit the observations by taking 
into account only the opacity provided by the Nai and 
K I resonance doublets and the continuum opacity sources 
listed in Table §. 
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Fig. 8. Gl 229B observed data (thin full line) compa red to 
some predicted spectra obtained using the Tsuji's ( 2000 ) 
C-type model of T c g ~ 1000 K and logg = 5. Different val- 
ues of a are oveplotted to show the effect of dust opac- 
ity at this low temperature. The best fit is provided by 
a = 0. 1 and N — 4. Only Cs I features are seen in the wave- 
lenth range presented here. 



Since the formation of dust in these cool atmospheres 
can produce additional opacity (AdO) which may affect 
the synthetic spectra we decided to investigate whether 
a simple description of it could help us to improve the 
comparison between observed and computed spectra. We 
adopted as law for AdO the following expression: a v = 
a Q * (v/vq) n . For A = to 4 this law corresponds to the 
case of radiation scattering produced by particles of differ- 
ent sizes, being N = 4 the case of pure Rayleigh scattering, 
and N = corresponding to the case of white scattering. 
However, at present we cannot distinguish whether this 
AdO is due to absorption or scattering processes. The pa- 
rameters N and do would be determined from the com- 
parison with observations, but in all cases we try to get 
the best fit for iV = 4, which would be the most simple 
from the physical point of view. We adopted as v Q the fre- 
quency of the K I resonance line at 769.9 nm. Our model of 
AdO is depth independent and therefore in this approach 
we cannot model the inhomogeneities (e.g. dust clouds) 
which may exist in L-dwarf atmospheres. 

We have investigated how our simple approach for the 
modelling of AdO may lead to better comparisons between 
predicted and observed spectra. The implementation of 
our law of AdO depresses the fluxes in the blue wing of 
the Ki doublet, considerably improving the reproduction 
of the observed data of Kelu 1 (Fig. |), DenisP J1228-1547 
(Fig. 0, upper panel) and DenisP J0205-1159 (Fig. @, lower 
panel). For each object, synthetic spectra were computed 



for a range of T e ff, gravities, ao and depletion factors R 
for TiO, VO, CaH and CrH. The previous figures display 
those syntheses which better reproduce the observed spec- 
tra. From our computations we infer that VO is less effi- 
ciently depleted than TiO at a given T e g. The depletion 
of these oxides appears to increase very rapidly as we go 
from BRI 0021-0214 and Kelul to the lower temperature 
objects DenisP J0205-1159 and DenisP J0205-1159. Note 
the good fit to the observed spectra at 860 nm provided 
by the (0,0) band of the A 6 £(+)-X 6 £(+) system of CrH. 
From the study of this band we also find that the depletion 
factor of CrH increases from the warmer to the coolest ob- 
jects in the sample. Finally, we also find that atomic lines 
become weaker and narrower when increasing the amount 
of AdO in the atmospheres. 

We have also studied whether we can explain the opti- 
cal spectrum of Gl 229B using the following Tsuji's model: 
C-type, T cff = 1000 K, and log g =5.0, and the AdO law of 
index N = 4 used above. In Fig. || we show several spectral 
synthesis reproducing the optical spectrum of this object, 
and showing the effect of different a parameters which 
is related to the amount of dust in the atmosphere. In 
G1229B we need the highest value of ao, which is inter- 
preted as evidence for the most "dusty" atmosphere in our 
sample. We have not attempted to reproduce the Cs I lines 
in G1229B. According to our hypothesis, the inclusion of 
AdO avoids the contribution of high pressure regions to 
the formation of these lines. The shorter wavelength Csi 
line at 894.3 nm is very affected by dust opacity. 



5. Discussion 

5.1. T c ff for L-type dwarfs 

Our computations provide a reasonable description of the 
far-red optical spectra of L-dwarfs and provide a physi- 
cal basis for a progressively decreasing T c g for th e pro- 
posed spectral classifications (Martin et al. 



1999; Kirk- 



patrick et al. 1999a ). Effective temperatures for a few field 
L-dwarfs and for Gl 229B have been deriv ed from spec- 
traat IR wavelength s (All ard et al. 1996 ; Marley et al. 



1996 ; Matthews et al. |1996j; Tsu ji et al. |1999| ; Jones et al. 
1996; Kirkpatrick et al. 1999b). Here we study to what 



extent we can use the broad energy spectral distribution 
in the optical to infer the T e g for the cool dwarfs in our 
sample. Our best estimates using Tsuji's models (see Ta- 
ble ^) are in good agreement with those found from IR 
data for objects of similar spectral types and for Gl 229B. 
Using Allard's models and the simple approach of AdO 



described in section 4.3 we are also able to reproduce the 



observed spectra, albeit we require lower T G ff's by up to 
several hundred degrees for the coolest L-dwarfs. This is 
mainly due to the hotter stratification of Allard's model 
atmospheres for the potassium and sodium lines forming 
layers (see Fig. ||). Opposite to IR-based temperature de- 
terminations, the estimation from optical spectra is very 
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Object 


Sp. Type 


I - J 


Teff 


Other measures 


Source 








(±200 K) 






BRI 0021-0214 


M9.5 


3.30 


2200 


1980, 2300 


TMR93, LAH98 


Kelul 


L2 (L2) 


3.50 


2000 


2000, 1900 


B99, RLA97 


DenisP J 1228-1547 


L4.5 (L5) 


3.81 


1600 


1800 


B99 


DenisP J0205-1159 


L5 (L7) 


3.82 


1200 


1700, 1800 


B99, TK99 



NOTES. Spectral types are given in Martin et al. (1999). Those spe ctral t ypes in brackets come from Kirkpatrick et al. (1999a 
(/ — J) colors have been taken from Leggett, Allard & Hauschildt ( 1998| ) 
Refere nces: TMR93 
et al. fll997| ); TK99 = 



: Tinney, Mould & Reid (|l993|); LAH98 = Leggett et al. ( 
Tokunaga & Kobayashi (|l999D . 



1998 



B99 = Basri et al. ( |l999| ); RLA97 = Ruiz 



sensitive to the input physics parameters, like dust forma- 
tion, molecular equilibrium, sources of opacity and atmo- 
spheric models, which limits the accuracy of the estimates 
to ^200 K. On the other hand, the optical spectra pro- 
vide an opportunity to test the reliability of the physical 
description of the atmospheres. 

Cs I lines in the optical spectra have been recently used 
to infer the T e s of some L-dwarfs (Basri et al. 1999). In 
spite of the sensitivity of these lines to many input param- 
eters in the models (i.e. to the amount of opacity, chemi- 
cal equilibrium, etc.) we find in general a good agreement 
with their estimated temperatures for the earlier L-dwarfs 
(see Table |]). However, for the latest L-type objects in 
our sample, DenisP J1228-1547 and DenisP J0205-1159, 
we find temperatures up to 400 K lower. We may attribute 
this discrepancy to the effects that dust opacity has in the 
formation of optical lines of alkalie. This produces a kind 
of "veiling" of the atomic lines, reducing their intensities. 
Basri et al. (|1999| ) did not consider this effect, and there- 
fore they required hotter models in order to explain the 
strength of the observed Cs lines. Another possible reason 
for the discrepancy is that the temperatures in Table [| 
have been obtained for gravity logc/ = 5; if we increased 
gravity by 0.5 dex we would have to increase the tempera- 
tures of our models by about 200 K. In this case the spec- 
tral synthesis does not reproduce the observations so well 
but they are still acceptable. 
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Fig. 9. Predicted optical spectra (Tsuji's C-type model 
atmosphere, T e ff= 1200K and logg = 5) for intermediate 
objects between G1229B and DenisP J0205-1159. Differ- 
ent intensities in the dust opacity have been considered. 
Spectra have been normalized at 860 nm. Identification of 
the atomic and molecular bands included in the computa- 
tions are also provided. 



5.2. Alkali lines: the case of lithium 

Recently, several new G1229B-like objects have been dis- 
covered by the SDSS survey (Strauss et al. 1999) and the 
2MASS collaboration (Burgasser et al. 1999). Based on 



their ncar-IR spectra, these authors suggest that these 
new cool brown dwarfs may be warmer than G1229B. In 
order to estimate the optical properties of these objects 
we have computed synthetical spectra using a Tsuji's C- 
type model of T c s = 1200K and logg = 5.0, and we have 
adopted the basic prescriptions that were followed in sec- 
tion 4, i.e total depletion of VO and TiO and the AdO law 
with N = 4. In Fig. || we plot the resulting spectra con- 
sidering different amounts of dust opacity. As expected, 



the overall shape of the spectrum is intermediate between 
that of DenisP J0205-1 159 and G1229B. In the absence 
of dust absorption (a o =0.0) alkali lines are clearly seen 
(including the lithium resonance doublet), and the spec- 
trum is governed by the sodium and potassium lines. If we 
consider dust opacities comparable to those in G1229B, 
the alkali lines of Cs and Rb become weaker, but still 
detectable with intermediate resolution spectroscopy. Re- 
markably, the subordinate Nal doublet at 819.5 nm is very 
sensitive to the incorporation of dust opacity due to the 
larger depths of its formation as compared with resonance 
lines. For very high dust opacities these lines may become 
undetectable. 
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Fig. 10. Densities of lithium and the most abundant 
species containing lithium are plotted for the C-type 
Tsuji's fl2000| ) model atmospheres with T cff = 2000K (up- 
per panel) and 1200 K (lower panel) and logp = 5. Com- 
putations have been performed under the assumption of 
complete chemical equilibrium for a solar elemental mix- 
ture. 



The effects of additional dust opacity on the formation 
of Li I lines (resonance and subordinate ones) also deserve 
detailed consideration since they play a major role as a 
discriminator of substellar nature for brown dwarf can- 



didates (see Rebolo, Martin, & Magazzu 1992; Magazzu, 
Martin, & Rebolo 1993 ). Most of the known brown dwarfs 
are actually recognized by the detection of the Li I reso- 
nance doublet in their spe ctra ( Rebolo e t al. 
et al. |l997a|; Rebolo et al. |l99S| ; Tinney |l998 ; 
et al. 1999a). The chemical equilibrium of lithium con- 



1996; Martin 



Kirkpatrick 
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Wavelength (nm) 

Fig. 11. Fitting of the Lil resonace line of Kelul 
(full line) using the Tsuji's C-type model atmosphere 
for T c ff = 2000K, logg = 5. Computations have been per- 
formed for a lithium abundace of log iV(Li) = 3.0 and con- 
sidering different amounts of dust opacity. The resolution 
of all theoretical spectra is the same than the one of the 
observed data. The predicted model that better fits the 
observations coincides with the one that also nicely repro- 
duces the overall shape of the optical spectrum (a Q = 0.03, 
see Fig. ||). 



tained molecules have been considered in all our synthe- 
ses (Fig. |Io| depicts the density profile s for lithium species 
for two C-type models by Tsuji 200C ). Our computations 
show that both, the resonance line at 670.8 nm, and the 
subordinate lines at 601.3 nm and 812.6 nm are very sen- 
sitive to the AdO that we need to incorporate in the spec- 
tral synthesis if we want to explain the observed broad 
spectral energy distribution. Among the subordinate lines 
the doublet at 812.6 nm is more easily detectable, but 
the predicted EWs, assuming fully preserved lithium, are 
rather small, ranging from EW = 0.4A to 0.04A for T e g- 
values in the range 2000 K down to 1200 K. These EWs 
are considerably reduced by the inclusion of the AdO 
described in the previous section, which makes their de- 
tection rather difficult. In Table || we give the predicted 
EWs of the Li I resonance doublet at 670.8 nm for several 
of the coolest model atmospheres (2000-1000 K) consid- 
ered in this work. First, we note that in the absence of 
any AdO (second column in the table), we would expect 
rather strong neutral Li resonance lines in the spectra of 
objects as cool as DenisP J0205-1159 and G1229B. The 
chemical equilibrum of Li-contained species still allow a 
sufficient number of Li atoms to produce a rather strong 
resonance feature; one reason for this is that CI and O 
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Table 5. Equivalent widths (A) of the Lil reso- 
nance doublet at 670.8 nm computed for the C-type 
Tsuji's ( 200C| ) model atmospheres, cosmic Li abundance 
(logiV(Li) = 3.2) and gravity log g = 5.0. 
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0.00 


0.01 


0.10 


(K) 




EW (A) 




1000 


17 


8 


0.6 


1200 


30 


12 


0.7 


1400 


42 


21 


0.9 


1600 


40 


24 


1.6 


2000 


23 


16 


3.6 



atoms should be also bounded into other molecules (e.g. 
NaCl, KC1, H 2 0, etc.). Our computations indicate that 
objects like DenisP J0205-1159 and cooler objects with 
moderate dust opacities should show the Li I resonance 
doublet if they had preserved this element from nuclear 
burning, and consequently the lithium test can still be ap- 
plied. Furthermore, even in very dusty cool atmospheres 
like that of Gl 229B for which we have inferred a high value 
of the opacity parameter of a = 0.1 (fourth column in Ta- 
ble ||), the lithium resonance line could be detected with 
an EW of several hundred mA (high S /N data would be 
required). 

Another effect that we shall consider is whether small 
changes in the AdO (which could be originated as a conse- 
quence of some "meteorological" phenomena occurring in 
these cool atmospheres) can lead to detectable variations 
in the EWs of the lithium lines. In particular, weak lithium 
lines do not necessarily imply a depletion of this element. 
The observed Li I variability in Kelu 1 (with changes in 
EW by a factor 5), could be an indication of meteoro- 
logical changes in the atmosphere of this rapidly rotating 
cool object. In Fig. [ll]we present several spectral synthesis 
showing the sensitivity of the lithium line to the AdO in 
the atmosphere. The AdO parameters which give the best 
fit for the lithium line in Kelu 1 (EW = 6.5±1.0 A) coincide 
with those also providing the best fit to the whole opti- 
cal spectrum (see Fig. ^). Anyway, the obtained lithium 
abundance is consistent with complete preservation of this 
element. 



6. Conclusions 

In this paper we have attempted to model the far-red spec- 
tra (640-930 nm) of several L-dwarfs suitably selected to 
cover this new spectral class. We have used model atmo- 
spheres from Tsuji ( |2000D and Allard ( |1999[ ), as we ll as an 
LTE spectral synthesis code (Pavlenko et al. 199E ) which 
takes into account chemical equilibrium for more than 100 



molecular species, and detailed opacities for the most rele- 
vant bands. We have arrived to the following conclusions: 

1) Alkali lines play a major role governing the far-red 
spectra of L-dwarfs. At early types, this role is shared with 
TiO and VO bands, which dominate this spectral region in 
late M-dwarfs. As we move to later spectral types we need 
to incorporate progressively higher depletions of these ox- 
ides and of the hydrides CrH and CaH, consistently with 
the expectation that Ti and V atoms are depleted into 
grains; and we also require additional opacity to reproduce 
the overall shape of the spectra. This additional opacity 
could be either due to molecular/dust absorption or to 
dust scattering. 

2) We have shown that a simple law for this addi- 
tional opacity of the form a {v/vo) , with N = 4, gives 
a sufficiently good fit to the observed spectra of L-dwarfs 
and G1229B. For this late object we require the highest 
value of a consistent with a very dusty atmosphere. The 
strength of alkali lines is highly affected by this opacity. 

3) From the best fits to our spectra, we derive the 
most likely T c g values for our sample of L-dwarfs. For 
the warmer objects, our values are consistent with those 
obtained by other authors, however we find lower T e ff's by 
serveral hundred degrees for the coolest L-dwarfs. Because 
the optical spectra are very much affected by the input 
physics, a more reliable T cf j scale should be obtained by 
fitting the IR data of these cool objects. 

4) After detailed consideration of chemical equilib- 
rium, we find that the lithium resonance doublet at 
670.8 nm can be detected in the whole spectral range 
(down to 1000 K). In the coolest L-dwarfs the strength 
of the resonance line is more affected by the amount of 
additional opacity needed to explain the spectra than 
by the depletion of neutral lithium atoms into molecular 
species. In those atmospheres where the additional opac- 
ity required is low, the lithium test can provide a useful 
discrimination of substellar nature. Changes in the phys- 
ical conditions governing dust formation in L-dwarfs, will 
cause variability of the lithium resonance doublet. Tak- 
ing into account the need for additional opacity in Kelu 1, 
we find that the lithium abundance can be as high as 
log iV(Li) = 3.0, i.e. consistent with complete preservation. 
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